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ABSTRACT 

In the near future there will be launched space missions (e.g. COROT, KEPLER), 
designed to detect Earth-like extrasolar planets. The orbital elements of these (still 
hypothetic) planets will contain some uncertainties, that can only be eliminated by 
careful dynamical investigations of the hosting planetary systems. The proportion of 
extrasolar planetary systems with one known giant planet is high (~ 90%), therefore 
as a first step we have investigated the possible existence of terrestrial planets in these 
systems. In this paper a development of a stability catalogue of the habitable zones 
of exoplanetary systems is reported. This catalogue is formed by a series of stability 
maps, which can help to establish, where Earth-like planets could exist in extrasolar 
planetary systems having one giant planet. After a description of the dynamical model 
and the numerical methods, details of the stability maps are discussed. An application 
of the stability catalogue to 15 known exoplanetary systems is also shown, and a 
characterization of the stability properties of their habitable zones is given. 
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1 INTRODUCTION 

Since the discovery of the firs t extrasolar plane t arou nd a 
, solar-type star, 51 Pegasi by iMavor fc Quelozl (|l995t ). by 
now (October 2006) 209 exoplanets have been found. These 
planets form 179 exoplanetary systems, among which there 
are 21 multiple systems. The main detection technique of ex- 
oplanets is based on the Doppler shift of the hosting star's 
spectral lines. This method is adequate to discover giant 
planets, but it is not yet capable to detect Earth-sized plan- 
ets. 

However, the detection and investigation of Earth- 
like exoplanets can be of great importance for the fur- 
ther improvements of formation theories of our Solar Sys- 
tem, where both rocky Earth-like, and gaseous Jupiter- 
like planets have been formed. Another major question 
which has a special interest not only for the scientific 
community is whether life may have evolved in some ex- 
oplanetary systems. The region of an exoplanetary sys- 
tem, where water-based life might evolve, is called the 
habitable zone (HZ). Simply saying, the HZ of a plane- 
tary system is that region , where liquid water can exist 
on th e surface of a planet (|Kasting. Whitmire fc Reynolds! 
Il993h . The dynamical stability of hypothetical Earth-like 
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planets in the HZs of several extrasolar planetary sys- 
tems has been the subject of intensive research (see for in- 
stance 
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Jones. Sleep fc Underwoodf 2006T ). According to these stud- 
ies many systems can host terrestrial exoplanets in their 
HZs. 

Recently, there have been much efforts for the construc- 
tion of space probes aiming at the detection of terrestrial ex- 
oplanets in the near future. Such is COROT, to be launched 
in November 2006, and KEPLER with a launch in 2008. 
These missions will use transit photometry as detection tech- 
nique. The transit of an unseen exoplanet in front of the 
hosting star's disc is expected to result in an observable dim- 
ming of the star's light intensity. There are also other plans, 
such as NASA's TPF and SIM, or ESA's Darwin and Gaia 
missions in the next decade aiming at studying all aspects 
of exoplanets. 

By using data obtained from the observations, the or- 
bital parameters of an exoplanet can be derived only within 
some error limits. In order to improve the orbital parame- 
ters certain constraints should be used. For instance, if an 
Earth-like planet were discovered in an extrasolar planetary 
system, where an already detected giant planet is known to 
orbit the star, such a constraint could be the dynamical sta- 
bility of the terrestrial planet. This means, that those orbital 
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solutions for the terrestrial planet, which result in unstable 
orbits, are unlikely and therefore should be excluded. 

The dynamical stability of Earth-like planets can be es- 
tablished in several ways, for example (i) by exploring the 
stable and unstable regions of the phase space of each ex- 
trasolar planetary system separately, or (ii) by using sta- 
bility maps computed in advance for a large set of orbital 
parameters. In this paper we present such stability maps 
and show how to apply them to the detected exoplanetary 
systems. The second method has the advantage that the 
stability properties of a terrestrial planet can be easily es- 
tablished when the orbital parameters of the giant planet 
of an exoplanetary system are modified. Very often the or- 
bital parameters of the giant planets are uncertain, and due 
to the accumulation and improvements of the observational 
data, they are changed quite frequently. If the first approach 
is chosen, one has to re-explore the phase space of the in- 
dividual exoplanetary system after each modification of the 
orbital parameters of the giant planet. However, this is not 
necessary when the second method is used, since the sta- 
bility properties of the investigated planetary system can 
be easily re-established from the already existing stability 
maps. 

As a first step towards general stability maps, we com- 
piled a stability catalogue for a dynamical model consisting 
of a star, a giant planet, and a small Earth-like planet. We 
present this catalogue in this paper. Section 2 describes the 
dynamical model and the applied numerical methods. De- 
tails of the catalogue and the stability maps are given in 
Sections 3 and 4. Section 5 shows how to apply the cata- 
logue to individual exoplanetary systems. By using this cat- 
alogue, the dynamical reality of the orbital elements of a 
newly discovered terrestrial planet in an extrasolar plane- 
tary system can be easily determined. This catalogue can 
also be used to establish the stability properties of the HZs 
of known exoplanetary systems. This is shown in Section 6. 



2 DYNAMICAL MODEL AND METHODS OF 
INVESTIGATION 

The majority of planetary systems, detected up to now, con- 
sists of a star and a giant planet revolving about it in an ec- 
centric orbit. Therefore we use a relatively simple dynamical 
model (known as the elliptic restricted three-body problem) , 
where two massive bodies, called the primaries, move in el- 
liptic orbits about their common center of mass, and a third 
body of negligible mas s moves under th eir gravitational in- 
fluence (for details see ISzebehely| [l967L In our application 
the primaries are a star and a giant planet, and the third 
body is a small Earth-like planet, regarded as massless. 

We note, however, that among the extrasolar planetary 
systems with long observational baselines, there is a high 
rate of multiple planet systems. Thus a more convenient 
model for our stability investigations would be the restricted 
iV-body problem with N - 2 giant planets (N > 4). The 
presence of additional giant planets may enhance the insta- 
bilities, induced by just one massive planet, and make the 
HZ of a system more unstable. The instabilities may appear 
due to mean motion resonances between the massive giant 
planet and the Earth- like planet. Thus, by mapping these 
resonances, we can find the most significant birthplaces of 



instabilities in the habitable zones. On the other hand, the 
dynamical model with one giant planet also offers the most 
convenient way to display the most important mean motion 
resonances as a function of the mass ratio of the star and 
the giant planet, and of the eccentricity of the giant planet. 
In this paper we mapped these resonances, and the effects 
of additional giant planets will be studied in a future work. 

To explore the stable and unstable regions of the mo- 
tion of the third body, the equations of motion have to be 
integrated for a large set of initial conditions and for a long 
time for each orbit. This can be very time consuming, even 
in the case of a huge computing potential. In order to study 
efficiently the stability properties of dynamical systems, sev- 
eral numerical techniques and chaos detection methods have 
been developed recently. Using the working hypothesis, that 
chaotic behaviour means instability through chaotic diffu- 
sion of trajectories in the phase space, these methods can 
be applied to explore the stable and unstable regions of 
dynamical systems. We are, however, aware of the fact, 
that chaotic behaviour does not necessarily mean instabil- 
ity. There could be orbits exhibiting chaotic behaviour, but 
stable in Nekhoroshev-sense, e.g. remaining in a bounded re- 
gion of the phase space for time intervals equal to the age of 
the investigated planetary system. However, if other desta- 
bilizing effects, such as the presence of yet undetected giant 
planets are present in the exoplanetary systems, chaotic be- 
haviour may result in really unstable orbits. In other words, 
in high dimensional dynamical systems unstable orbits em- 
anate only from the chaotic regions of the phase space. Ex- 
ploring these regions we can find the birth domains of unsta- 
ble orbits. However, in order to decide the stability character 
of an orbit, one has to use long-term numerical integration of 
the equations of motion. On the other hand, in the case of a 
large set of initial conditions the fast chaos detection meth- 
ods allow us to explore qualitatively the dynamical structure 
of the extrasolar planetary systems. 

Our investigations are based on such chaos 
detection methods. As a main tool we used the 
method of the relative Lya pu nov i ndi cators (RLI) 
llSandor. Erdi fc Efthvmiopoulosl l200d ; ISandor et al.1 
l2004h. and in some cases the fast Lyapunov indicators (FLI) 
llFroeschle. Lega fc Gonczil 1 19971 ; iGuzzo. Lega fc Froeschli 
2002). We also applied t he maximum eccent ricity method 
(MEM) (see for instance iDvorak et al.ll2003T l. as an inde- 
pendent check of the RLI and the FLI by direct numerical 
integration of the equations of motion. The knowledge of 
the eccentricity of stable orbits is essential for the problem 
of the development of a biosphere on a terrestrial planet in 
the habitable zone of a star. For example, recent estimations 
for the maximum eccentricity for the Earth - after detailed 
atmosphere modelfl- are e ma x < 0.2. 

Regarding the RLI and the FLI, our choice is motivated 
by the complementary character of the two chaos detection 
methods. The RLI measures the difference between the con- 
vergence of the finite time Lyapunov indicators to the max- 
imal Lyapunov characteristic exponent of two initially very 
close orbits. This method is extremely fast in determining 
the ordered or chaotic nature of individual orbits, even to 
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ditions the RLI (or FLI) value of the corresponding or- 
bit of the test planet is computed and visualized in the 
(a, ei) parameter plane. Due to the very small stepsize in 
a and ei, each stability map corresponds to more than 
8 x 10 4 initial conditions, thus providing a very fine res- 
olution. In the following section we show some stability 
maps from the catalogue, containing the HZs of some ex- 
oplanetary systems as examples. The complete catalogue 
is available at http://astro.elte.hu/exocatalogue and 
http : //www. astro .univie . ac . at/ adg. 



4 DETAILS OF THE STABILITY MAPS 

In Fig. [TJ we displayed a case, where the semi-major axis 
of the test planet's orbit is larger than that of the giant 
planet (a > 1, and /i = 0.001, Mi = 0°). Light regions corre- 
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detect orbits with very s mall Lyapunov cha racteristic ex- 
ponents (e.g. sticky orbits. ISandor et al.ll2004h . as well as to 
distinguish between ordered and chaotic regions of the phase 
space. According to our experiments, gained in different dy- 
namical problems, it is enough to integrate the two initially 
close orbits for a few hundred times of the longest orbital 
period of the investigated system. As to the FLI, their time 
evolution shows clearly the ordered or chaotic character of 
an orbit, and the method is very effective in detecting sta- 
ble resonant orbits. Concerning the computation of the FLI, 
the integration time is several tens of thousands orbital peri- 
ods. In the MEM we compute by numerical integration the 
largest value of the orbital eccentricity, reached by an or- 
bit during a given, long time interval (10 5 orbital periods in 
this paper). These methods have been successfully applied 
to several exoplanetary systems to determ ine their dynam- 
ical stability prope rties: for the RLI see (lErdi et al.l 12 004; 

Sandor et all 12004). for the FLI dPilat-Lohinger fc Dvorak! spom } to low values of the RLI {W~ 10 ), and thus to ordered, 

dynamically stable motion of the test planet. Dark shades 
correspond to large values of the RLI (10 -5 ) and chaotic be- 
haviour of the test planet. The dark unstable regions on the 
left are due to the proximity of the giant planet. The fig- 
ure is dominated by V-shaped gray stripes, corresponding 
to outer mean motion resonances between the test planet 
and the giant planet. These resonances, marked at the top 
of the figure, can represent either ordered (stable for infinite 
time), or weakly chaotic (which may become unstable after 
very long time) behaviour. With the increase of e\ many res- 
onances overlap, giving rise to strongly chaotic and thus very 
unstable behaviour. The reason for this is that by increasing 
ei the apocenter distance of the giant planet also increases, 
and the giant planet perturbs more strongly the outer test 
planet. The figure also shows, for two planetary systems, the 
positions of the zero age HZ, which we call the classical HZ, 
and the present HZ. We recall, that the habitable zone (HZ) 
is that region around a star, where water can exist in liquid 
phase on the surface of a planet by virtue of stellar irradi- 
ance. The boundaries of the classical HZ are calculated by 
assuming the physical properties of a zero age main sequence 
star, while the boundaries o f the present HZ d epend also on 
the age of the hosting star (I Jones et al.ll2006t) . The borders 
of the classical HZs for a large number of extrasolar plane- 
tary systems are given bv lMenou fc Tabachnikl (|2003l ). while 
the borders of the presen t HZs can be found in the more re- 
cent work of Ijones et al.1 (|2006t ). We note that according to 
the latter work the present HZs are shifted outwards with 
respect to the location of the classical HZs. 

Fig. [5] shows a stability map for the same initial condi- 
tions as in Fig. [TJ but resulting from the FLI computations. 
It has the same overall structure, as the map obtained by 
the RLI. Note, however, the complementary character of the 
results obtained by the two methods. The RLI detects the 
chaotic separatrices of the resonances (dark V-shaped stripes 
in Fig. [TJ, while the FLI finds the stable resonant orbits in- 
side the resonances (light stripes inside the resonances in 
Fig. [2}. Note, that for the given initial conditions (Mi = 0°) 
almost all resonances can have stable orbits, but the 1:3 res- 
onance is chaotic according to bot h methods, a w ell known 
result for the 3:1 inner resonance (|Wisdomlll983h . 

Fig. |3] has been calculated by the MEM, again for the 
same initial conditions as in Figs. [TJ and [2] It can be seen, 
that similarly to the RLI and the FLI, the MEM also indi- 



3 CATALOGUE OF DYNAMICAL STABILITY 

We compiled a catalogue of dynamical stability, adapted to 
a system consisting of a star with mass mo, a giant planet 
with mass mi, and a small test planet with negligible mass. 
The total mass of the star and the giant planet is taken as 
unit mass, thus /i = mi/(mo + mi) is the mass parameter of 
the problem, which has been changed between (i) 1 x 10~ 4 
and 9 x 10~ 4 in steps of A/i = 10 -4 , (ii) between 1 x 10~ 3 
and 9 x 10 -3 in steps of A/z = 10 -3 , and (iii) between 1 x 
10~ 2 and 5 x 10~ 2 in steps of A/i = 10~ 2 . This results 
in 23 values for fi. The giant planet was assumed to move 
around the star in an elliptic orbit, whose plane is taken 
as reference plane, with semi-major axis oi, eccentricity e%, 
argument of periastron ui, and mean anomaly Mi. Without 
the loss of generality, the semi-major axis a\ can be taken 
as unit distance, thus oi = 1. We assumed that wi =0°. 
The eccentricity ei was changed between 0.0 and 0.5 with 
a stepsize of 5 x 10~ 3 . During our computations we also 
changed Mi between 0° and 360° with A Mi = 45° both 
for inner and outer orbits of the test planet (with respect to 
the orbit of the giant planet). Additionally, we changed Mi 
for inner orbits of the test planet between 0° and 90° with 
AMi = 10° too. We assumed that the test planet moves in 
the orbital plane of the giant planet, and we always took 
e = 0, uj = 0°, M = 0° as initial orbital elements of the 
test planet. The case when the orbit of the test planet is 
inclined with respect to the giant planet's orbit will also be 
studied in future work. We changed the semi-major axis a 
of the test planet in two different intervals: (a) for orbits 
of the test planet "inside" the orbit of the giant planet we 
changed a between 0.1 and 0.9 with a stepsize of 10 -3 , and 
(b) for "outside" orbits between 1.1 and 4.0 with a stepsize 
of 3.625 x 10" 3 . 

The catalogue consists of 552 stability maps, corre- 
sponding to each (/j,, Mi) pair of the 23 values of fj, and 
16 values of Mi for inner orbits (368 maps) and 8 val- 
ues of Mi for outer orbits (184 maps) of the test planet. 
For a given (fj,, Mi) pair the stability map is prepared in 
such a way, that for each set of the described initial con- 
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cates the mean motion resonances. Since in the MEM orbits 
were integrated for 10 5 periods of the giant planet, we used 
it as a test of the two other chaos detection methods. 

Figs. U and [5] show stability maps for outer orbits of the 
test planet, computed by the RLI for /j, — 0.002 and /i = 
0.007, respectively. It can be seen, that for larger values of /j, 
the chaotic regions are more extended, due to the stronger 
perturbing effect of the giant planet. 

In Figs. [6] - 1101 we display cases where the semi-major 
axis of the test planet is smaller than that of the giant 
planet (a < 1). Two panels for fi = 0.001 (Figs. [6K and 
[6p) show the results for two different starting positions of 
the giant planet, Mi = 0° (Figs, (gja,) and Mi = 180° (Figs. 

). Here again mean motion resonances, this time inner 
ones, between the test planet and the giant planet domi- 
nate the stability maps. Inside the resonances the stable or 
chaotic behaviour of the test planet depends on the initial 
angular positions of the two planets. This is clearly visi- 
ble, when we compare the two panels. In Fig. [6Jl the giant 
planet starts from the pericenter, in Fig. |Sp from the apoc- 
enter, while the initial position of the test planet in both 
cases is at M = 0°, corresponding to the direction of the 
giant planet's pericenter. Comparing the two figures it can 
be seen that the location of the resonances does not change, 
since this depends on the ratio of the semi-major axes of 
the two planets. However, the character of the resonances 
becomes different. In Fig. [6p several resonances (5:2, 5:3, 
3:2) are more stable than in Fig. [6^,. This is due to the fact 
that the relative initial positions determine the places of con- 
junctions of the two planets. If they meet regularly near the 
pericentre of the giant planet, the motion of the test planet 
becomes chaotic, while it can remain regular if the conjunc- 
tions take place near the apocenter of the giant planet. The 
effect of the initial phase difference between the planets is 
important, therefore we prepared stability maps for more 
initial values of the mean anomaly Mi of the giant planet. 
The essential role of the phases for computations of orbits in 
mean motion resonances for terrestrial plane ts in exoplan- 
etary systems has already been stressed by lAsghari et al.l 
(2004), where the giant planet was initially in periastron or 
apastron, and the phase of the terrestrial planet has been 
set to 12 different initial angles of the mean anomaly. (A 
detailed treatment of the geometry of the resonances can be 
found in lMurrav fc Dermottll 19991 ). 

By studying the series of Figs. [S] - 1101 one can follow 
the effect of the increase of the mass parameter on the sta- 
bility structure. For larger mass of the giant planet, thus for 
larger perturbations, the chaotic region is more extended, as 
expected. 

We note also a peculiar feature of Figs.[6k,,[7]-[10] In the 
strongly chaotic regions on the right, which are due to the 
proximity of the giant planet, there are two almost parallel 
light stripes, corresponding to satellite-type orbits around 
the giant planet. These orbits are less chaotic (or they may 
be ordered), than those emanating from the surrounding re- 
gion. Between these stripes the dark gap corresponds to the 
perihelion distance of the giant planet, therefore in orbits 
originating from here, the test particles suffer close encoun- 
ters with the giant planet, or they even can collide with it. 



5 APPLICATION OF THE CATALOGUE TO 
INDIVIDUAL EXOPLANETARY SYSTEMS 

In this section we show how to use the catalogue to de- 
termine the stability of hypothetical Earth-like planets in 
exoplanetary sytems. As an example we consider the case 
of HD10697, where ai = 2.13 AU, ei = 0.11, ^ = 0.0055. 
Fig. [9] shows a stability map, calculated for /i = 0.005 for 
inner orbits of the test planet. This corresponds to the 
minimum mass of the giant planet (minimum masses are 
used throughout). The stability of a small planet (start- 
ing with e = 0) in the system HD 10697 can be studied 
along the line ei = 0.11. One can see, that for small semi- 
major axes, a < 0.33ai(= 0.729AU) the parameter space is 
very stable. When a > 0.33ai several resonances appear, 
among which the most important are the 5 : 1, 4 : 1, 
3 : 1, and 2 : 1 resonances. For a > 0.73ai(= 1.55AU) a 
strongly chaotic region appears. The classical HZ of this sys - 
tem is between 0.85 - 1.65 AU l|Menou fe Tabachnik|[2003T ). 
Thus in Fig. [9] the scaled classical HZ is situated between 
0.85/ai = 0.39 and 1.65/ai = 0.77 (shown as a rectan- 
gle, elongated in horizontal direction). One can see, that the 
inner part of the classical HZ contains ordered regions, but 
stripes of certain resonances are also present. The outer part 
of the classical HZ is in the strongly chaotic region. 

In Table [1] we list the orbital parameters of the ex- 
oplanetary systems displayed in Figs. [Tl HOI and we char- 
acterize the stability properties of their hypothetical ter- 
restrial planets in their cla ssica l and presen t HZs given by 
iMenou fc Tabachnikl <|2003h and l Jones et~aH (|2006l ). We give 
also the boundaries of the stable and the strongly chaotic 
regions. As sources for the parameters we used the web- 
pages http://www.exoplanets.org and Jean Schneider's 
catalogue (http : / /exoplanet . eu) , denoted by E and JS, re- 
spectively, in Table 1 and in the figures, when there were sig- 
nificant differences between the data. In the column MMRs 
we display the highest and the lowest order mean motion res- 
onances of the given exoplanetary system (not only in their 
HZ), the other resonances between them can be taken from 
the corresponding figure. When the highest order resonance 
is above 9 : 1, we call it as h.o.r. In the column CHZ/PHZ 
we give a characterization of the stability properties of the 
classical and present HZs. The notation S means a fully sta- 
ble HZ, where no mean motion resonance exists. When the 
HZ contains a few resonances, we consider it partly stable 
(PS). If many resonances fill in the HZ, we call it marginally 
stable (MS). Strongly chaotic, thus very unstable domains 
are denoted by SC. 

The uncertainty in the eccentricity ei of a giant planet 
can affect the stability of the HZ significantly. The uncer- 
tainty in ei is indicated on the stability maps by the width of 
the rectangle of the classical HZ. (In the case of the present 
HZ, marked by parentheses, we did not show the uncertainty 
in ei, since it is the same.) When the uncertainty in e\ is 
large, the stability properties at the top and at the bottom 
of the rectangle of the classical HZ can be quite different. 
For instance, in the case of the system HD 30177 it can be 
seen from Fig. 10, that for small values of e\ there are a few 
resonances in the classical HZ, thus it is partly stable, while 
for high values of ei the HZ is full of resonances, therefore 
it is only marginally stable. 



A stability catalogue of the habitable zones 5 



6 COMPARISON WITH OTHER RESULTS 



The habitability of individual exoplanetary systems has 
been studied by many authors. As we have already men- 
tioned, the orbital parameters of these systems are changed 
quite frequently, so it is not an easy task to compare di- 
rectly our results with the results of other authors using 
older or different orbital data. On the other hand, there are 
only a few works studying the habitability of a large portion 
of the k nown systems. Such g l obal i nvestigations have been 
made bv lMenou fc Tabachnikl (|2003l ) studying the habitabil- 
ity of 85 kn o wn ex oplanetary sytems, and more recently by 
IJones et al.l l|2006l ) investigating the present habitability of 
152 exoplanetary sytems. 

In these studies the gravitational influence of the known 
giant planet on the HZ is estimated mainly by some mul- 
tiple of its Hill radius. However, the stability properties of 
the HZ of an exoplanetary system are also affected by mean 
motion resonances. If the HZ of a system contains only a few 
resonances, the ordered regions between them can provide 
stability for an Earth-like planet. If the resonances cover 
densely the HZ, the long term stability of an Earth-like 
planet is more questionable; either the resonances overlap 
each other generating a strongly chaotic region, or chaotic 
diffusion, acting inside a mean motion resonance, can desta- 
bilize the motion of a hypothetical small planet. We note 
that the presence of additional, yet undetected giant plan- 
ets in an extrasolar planetary system can enhance the deple- 
tion of the mean motion resonances through their possible 
secular interactions with the known giant planet, similarly 
to the formation of the Kirkwood gaps in our Solar s ystem 
l|Morbidelli fc Moons! [l993l ; iMoons fc Morbidelli|[l995l ). 

In our stability maps we investi gated the dynamical 
stabili ty of the classical HZ given by iMenou fc Tabachnikl 
|2003t ) , and the present HZ calculated by Jones et al.l 1 2006(h 
Studying the dynamical stability of the classical and the 
present HZs of 15 extrasolar planetary systems (displayed in 
Table 1 ), we have found quite good agr eement with the re - 
sults of IMenou fc Tabachnikl (120031 ) and I Jones et al.l (|2006t) . 
We emphasize, however, that our stability maps also exhibit 
the resonant structure of the HZs, which is essential to es- 
tablish the real stability properties of an Earth-like planet. 

In what follows we characterize the HZ of some individ- 
ual exoplanetary systems. 

HD 52256 and HD 121504: Depending on the data source, 
the classical HZ of the system HD 121504 is either fully sta- 
ble (source E), or partly stable (source JS), containing a few 
resonances (Fig.[T|. The classical HZ of the system HD 52256 
is partly stable, but contains more resonances and smaller 
stability regions for the data JS, than for E. The present 
HZs of these systems are located in the very regular regions 
of the (a,ei) parameter plane with a few resonances in the 
case of HD 52256 for the data JS. Our results are in good 
agreement with the results of M cnou fc Tabachnikl (|2003l ) 
and I Jones et al.l (|2006l ) concerning the extension of the sta- 
bility regions. These systems could harbor stable habitable 
Earth-like planets. 

HD 8574: The classical HZ of this system is in the strongly 
chaotic region for both data sources E and JS (Fig. [4]). 
IJones et all l|2006h found that 72% of the present HZ of the 
system can survive the gravitational influence of the giant 
planet, thus it could contain a stable habitable planet. Our 
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Figure 1. Stability map, computed by the RLI for outer orbits of 
the test planet, when fi = 0.001 and Mi = 0°. The classical HZs 
of two exoplanetary sytems (with orbital parameters taken from 
two different data sources, denoted by E and JS, see Section 5) are 
also indicated as elongated rectangles, whose width corresponds 
to the uncertainty in the eccentricity of the giant planet's orbit of 
each system. The inner borders of the present HZs are displayed as 
opened parentheses. The outer borders are out of the investigated 
range of the (a, ei) plane. As described in Sections 5 and 6, it can 
be seen, that depending on the data source the classical HZ of 
the system HD 121504 is either fully stable (source E), or partly 
stable (source JS), containing a few resonances. The system HD 
52265 is partly stable, but contains more resonances and smaller 
regions of stability for the data JS, than for E. The present HZs 
are located in the very regular regions of the (a, ei) plane. 
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Figure 2. Stability map, computed by the FLI for outer orbits of 
the test planet, when fi = 0.001 and Mi = 0°. Light shades indi- 
cate stable regions, dark means chaoticity. The very light stripes 
correspond to stable resonant orbits in the mean motion reso- 
nances. Note the complementary character of Figs. [T] and [2] de- 
scribed in the text. 



approach clearly shows that the present HZ of the system 
is full with overlaping mean motion resonances, and only 
partly stable for the data E, and marginally stable for the 
data JS. Thus the chance for the system to host a stable 
habitable planet is small. 

70 Virginis and HD 178911 B: The classical HZ of 70 Vir 
is partly in the strongly chaotic region (Fig.0, while its re- 
maining part is full with overlaping mean motion resonances 
and thus only marginally stable. This means that the classi- 



6 Sdndor et al. 



see http://astro.elte.hu/exocatalogue 
or 

http://www.astro.univie.ac.at/adg 



see http://astro.elte.hu/exocatalogue 
or 

http://www.astro.univie.ac.at/adg 



Figure 3. Stability map, computed by the MEM for outer orbits 
of the test planet, when /i = 0.001 and Mi = 0°. Here again 
lighter shades mean orbits, where the eccentricity remains rel- 
atively small. The changes in the maximum eccentricity of the 
test planet indicate clearly the positions of the mean motion res- 
onances. 



Figure 5. Stability map, computed by the RLI for outer orbits 
of the test planet, when fi = 0.007 and Mi = 0°. As described 
in Section 6, the classical HZ of 70 Vir is partly in the strongly 
chaotic region, while its other part is marginally stable. The clas- 
sical HZ of the system HD 178911 B is partly stable. The present 
HZ of HD 178911 B (denoted by parenthesis) is stable, while at 
70 Vir it is partly stable, containing higher order mean motion 
resonances (1:8, 1:9, etc., not shown in the figure). 
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Figure 4. Stability map, computed by the RLI for outer orbits of 
the test planet, when fi = 0.002 and Mi = 0°. A comparison with 
Fig. 5 reveals the increase of the chaotic region for increasing /i. 
As described in Section 6, the classical HZ of the system HD 8574 
is strongly chaotic for both data sources E and JS. The present 
HZ of the system (denoted by parentheses) is partly stable for E, 
and marginally stable for JS. 



cal HZ of this system is unlikely to host habitable Earth-like 
planets. The classical HZ of HD 178911 B is partly stable, 
and it can contain stable Earth-like pl anets. Our results are 
in goo d agreement with the results of iMenou fc Tabachnikl 
(2003). The present HZs of these systems are in the or- 
dered regions of the (a, ei) plane, thus these systems can 
host Earth-like habitable planets i n stable orbi t s. Th ese re- 
sults agree well with the results of iJones et al.l (|2006h - 

Eps Eri, HD 114729, and GJ 777 A: The classical HZs of 
Eps Eri and HD 114729 are marginally stable (Fig. |6j). 
There is only a small chance to find Earth-like planets here. 



IJones et all (l2006t ) found that only 17% of the present HZ 
of Eps Eri can have stable orbits, and nowhere in HD 114729 
does so. In the case of GJ 777 A, due to the large uncer- 
tainty in ei of the giant planet (ei = 0.48 ± 0.2), the classi- 
cal HZ can be stable (for ei = 0.48 — 0.2), or partly stable 
(for e\ = 0.48). There are no computed stability data for 
ei > 0.5, however, the stable region will certainly decrease 
for larger values of e\. The present HZ of GJ 777 A for 
ei = 0.48 is partly in the strongly chaotic region of the 
(a,ei) plane, while its remaining part is filled by overlap- 
ing mean motion resonances and is marginally stable. For 
ei = 0.28, t he pr esent HZ of GJ 777 A is partly stable. 
IJones et alj (|2006l ) found that 64% of the present HZ can 
survive the gravitational influence of the giant planet. How- 
ever, even though resonances were included in the develop- 
ment of their fast method of assessing stability (jjones et al.l 
2005), the method does not permit detailed consideration 
of the resonances in the individual systems to which it 
was applied. Considering the classical HZ, this system can 
host stable Earth-like plan ets, as has also been fou nd by 
IMenou fc Tabachnikl i|2005t ) and lAsghari et al.l (|2004T ). 

HD 70642: The classical HZ of this system can be stable, 
or partly stable, depending on the eccentricity of the giant 
planet (Fig. [7J). The present HZ of the system contains a 
few mean motion resonances, but between them there are 
ordered regions, which can host s table Earth- l ike p lanets. 
This result has also been found by I Jones et al. (2006). 

HD 72659: The classical HZ of this system contains only 
a few mean motion resonances (Fig. [8}, thus it can host 
stable Earth -like planets. This is in go od a greement with 
the re s ults oflMenou fc Tabachnikl (|200&t ) and lAsghari et all 
|2004h . |jones et al.li|2006r i found, that 71% of the present HZ 
is not disrupted gravitationally by the giant planet. This is 
true, if one considers only the effect of the strongly chaotic 
region. However, the inner part of the present HZ is full with 
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Figure 6. Stability maps, computed by the RLI for inner orbits 
of the test planet, (a) /j, = 0.001, Mi = 0°; (b) fj, = 0.001, Mi = 
180°. The comparison of the two figures shows the effect of the 
relative starting positions of the two planets (for more details see 
the text in Section 4). As described in Section 6, the classical 
HZs of Eps Eri and HD 114729 are marginally stable. Due to the 
large uncertainty of the giant planet's eccentricity, the classical 
HZ of GJ 777 A can be stable (for ei = 0.28), or partly stable 
(for ei = 0.48). The HZ of the Solar System is placed on the 
map under the approximate condition that there is only one giant 
planet, Jupiter in the system. It is reassuring that the HZ of the 
Solar System is in the stable region. The present HZ of GJ 777 A 
for ei = 0.48 is partly in the chaotic region, while its other part 
is marginally stable, and for e\ = 0.28 it is partly stable. 



mean motion resonances and is marginally stable. Thus the 
place for a stable Earth-like planet is much smaller, than 
71% of the whole present HZ, when the effects of the mean 
motion resonances are taken into consideration. However, 
the presence of an Earth-like planet cannot be excluded en- 
tirely. 

HD 50554: Both the classical and the present HZ of this 
system are entirely in the strongly chaotic region of the 
(a, ei) plane (Fig. [9j> . Thus this sy s tem c an not host habit- 
able Earth-like planet. I Jones et all l|200rj ) came to the same 
conclusion. On the other hand, there could be (retrograde) 
satellite-type orbits in the system. 
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Figure 7. Stability map, computed by the RLI for inner orbits 
of the test planet, when /i = 0.002 and Mi = 0°. Due to the large 
uncertainty of the eccentricity of the giant planet, the classical 
HZ of the system HD 70642 is either stable, or partly stable. The 
present HZ of this system is partly stable within the lower and 
upper limits of ei. 
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Figure 8. Stability map, computed by the RLI for inner orbits 
of the test planet, when /i = 0.003 and M\ = 0°. The classical 
HZ of the system HD 72659 is partly stable. The outer part of 
the present HZ of the system is chaotic, while its inner part is 
marginally stable. We placed also the classical HZ of the multiple 
planet system 47 UMa on the figure for comparison, by taking into 
account from its two giant planets the more massive component 
b. 



14 Her: Considering the classical HZ, this system can host 
stable Earth-like planets, however, mean motion resonances 
fill densely the classical HZ (Fig. [9}. As to the present HZ, 
the place for a stable Earth-like planet is much smaller, with 
one part in the chaotic region, and the other is filled by 
many resona nces. Our re s ults a re in good agreement with 
the results of ljones et"ai"1 (|2006r i. 

HD 10697: The classical HZ of this system contains a 
strongly chaotic region, while its inner part between the 
mean motion resonances contains ordered, thus stable re- 
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Figure 9. Stability map, computed by the RLI for inner orbits 
of the test planet, when fj, = 0.005 and Mi = 0°. Along the line 
ei = 0.11, corresponding to the system HD 10697, the parameter 
space is very stable for a < 0.33. For a > 0.73 the parameter 
space is strongly chaotic. Between 0.33 < a < 0.73 there are 
several resonances. The scaled classical HZ of the system HD 
10697 is between 0.39 < a < 0.77, thus its inner part is partly 
stable, but its outer edge is strongly chaotic. So is its present 
HZ. As described in Section 6, both the classical and the present 
HZ of the system HD 50554 are strongly chaotic. Interestingly, 
there could be (retrograde) satellite-type orbits in this system. 
The classical HZ of 14 Her is partly stable, while its present HZ 
is chaotic on one side, and marginally stable on the other. 



Figure 10. Stability map, computed by the RLI for inner orbits 
of the test planet, when /i = 0.009 and M\ = 0°. Due to the large 
uncertainty in the eccentricity of the giant planet, the classical HZ 
of the system HD 30177 is either partly stable (at the bottom of 
the rectangle), or marginally stable (at the top of the rectangle). 
The present HZ is strongly chaotic for large values of ex, while 
its inner part is partly or marginally stable. 



astronomers in determining the orbital parameters of ter- 
restrial planets detected in the future by space instruments 
(COROT, KEPLER) using transit photometry. 



gions as well (Fig. [9]l . The present HZ of this system is in 
the strongly chaotic region, thus it ca nnot harbour stabl e 
Earth-like planets. This is also found bv I Jones" et al] l|200fj ). 

HD 30177: The classical HZ of this system may host Earth- 
like planets, since between the mean motion resonances 
there are ordered regions on the (a, ex) parameter plane 
(Fig. llO[) . Considering the present HZ, its outer part is in the 
strongly chaotic region, while the inner part of it may con- 
tain stable Earth- like pl a nets. This result agrees well with 
the result of Ijones et al] l|2006l ). 

Thus we can conclude, that the effect of the mean mo- 
tion resonances cannot be neglected, when studying the hab- 
itability of the known extrasolar planetary systems. In this 
sense o ur work compleme nts the already existing studies, 
such as IJones et al] feoOfj ). 



7 SUMMARY 

In this paper we described a catalogue of stability maps 
of hypothetical terrestrial planets in exoplanetary systems 
with one giant planet. By using these maps, one can easily 
study the stability properties of possible Earth-like planets 
and the habitable zones in known or future exoplanetary 
systems. The novelty of our work is the invariance of our re- 
sults with respect to the frequent changes in the computed 
orbital parameters of the extrasolar planetary systems. Thus 
any new orbital determination (leading to new elements) can 
immediately be checked without any time-consuming com- 
putations. This catalogue could also be useful for observing 
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Table 1. Parameters of the extrasolar planetary systems displayed in Figs, H HTOl The semi-major 
axis a\ and the boundaries of the stable and the strongly chaotic regions are in AU. MMR means 
mean motion resonance, CHZ stands for the classical, PHZ for the present habitable zone. The 
characterization is: S stable, PS partly stable, MS marginally stable, SC strongly chaotic (n.d. 
means, that there are no data for the present HZ). The data sources E and JS are explained in 
the text. 



System 


ai 


ei 


stable 
region 


MMRs 


strongly chaotic 
region 


CHZ/ 
PHZ 


HD 52265 (E) 


0.49 


0.19 


a > 1.47 


2:3- 


+1:5 


a < 0.735 


PS/S 


HD 52265 (JS) 


0.49 


0.29 


a > 1.47 


2:3- 


+1:6 


a < 0.735 


PS/PS 


HD 121504 (E) 


0.33 


0.03 


a > 0.69 


2:3- 


+1:3 


a < 0.429 


s/s 


HD 121504 (JS) 


0.33 


0.13 
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2:3- 


+1:5 


a < 0.49 


PS/S 


HD 8574 (E) 


0.77 


0.288 
± 0.05 


a > 2.618 


2:3- 


+1:6 


a < 1.54 


SC/PS 


HD 8574 (JS) 


0.77 


0.4 
± 0.04 


a > 2.85 


2:3- 


+1:7 


a < 1.54 


SC/MS 


70 Vir 


0.48 
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2:3- 


+1:7 


a < 1.04 


MS/PS 


HD 178911 B 
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1:3- 


+1:6 
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PS/S 


Eps Eri 


3.4 


0.43 


a < 0.58 


h.o.r. 


-» 4:1 


a > 1.36 


MS/n.d. 


HD 114729 


2.08 


0.31 


a < 0.47 


8:1- 


+3:1 


a > 1.04 


MS/n.d. 


GJ 777 A 


4.8 


0.48 
0.48 - 0.2 


a < 0.82 
a < 1.10 


h.o.r-^7:3 
9 : 1,8 : 1 


a > 1.63 
a > 2.69 


PS/MS 
S/PS 


HD 70642 


3.3 


0.1+0.06 
0.1 - 0.06 


a < 0.89 
a < 1.32 


6:1- 
5:1, 


+2:1 
4:1 


a > 2.24 
a > 2.57 


PS/PS 
S/PS 


HD 72659 


4.16 


0.2 


a < 1.04 


9:1- 


+7:3 


a > 3.12 


PS/MS 


47UMa b 


2.09 


0.06 


a < 0.71 


5:1- 


+2:1 


a > 1.52 


PS/ 


HD 50554 


2.38 


0.42 
±0.03 


a < 0.24 


h.o.r- 


->5:1 


a > 0.83 


SC/SC 


14 Her 


2.85 


0.38 


a < 0.59 


h.o.r. 


-+4:1 


a > 0.99 


PS/MS 


HD 10697 


2.13 


0.11 


a < 0.72 


5:1- 


+2:1 


a > 1.4 


PS/SC 


HD 30177 


3.86 


0.3 + 0.17 
0.3 - 0.17 


a < 0.19 
a < 1.16 


h.o.r. 
h.o.r. 


-►7:1 
-+7:3 


a > 1.16 

a > 2.9 


MS/SC 
PS/PS 
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